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Noncollinear magnetic states in clusters are studied by using the single-band Hubbard Hamil- 
tonian. The unrestricted Hartree-Fock (UHF) approximation is considered without imposing any 
symmetry constraints neither to the size or orientation of the local magnetic moments {Si) nor to the 
local charge densities (ni). A variety of qualitatively different selfconsistent solutions is obtained as 
a function of cluster size, structure, number of valence electrons v and Coulomb interaction strength 
U It. This includes inhomogeneous density distributions, paramagnetic solutions, magnetic solutions 
with collinear moments and noncollinear spin arrangements that show complex antiferromagnetic 
and ferromagnetic-like orders. The environment dependence of the magnetic properties is analyzed 
giving emphasis to the effects of antiferromagnetic frustrations in compact structures close to half- 
band filling. Electron correlation effects are quantified by comparing UHF and exact results for the 
local magnetic moments, total spin, spin-correlation functions and structural stability of 13-atom 
clusters. Goals and limitations of the present noncollinear approach are discussed. 



PACS numbers: 36.40. Cg, 71.24.+q, 75.10.Lp 



I. INTRODUCTION 

Clusters may show specific phenomena that do not 
have an equivalent in the thermodynamic limit and that 
are therefore of interest in both basic and applied science. 
Moreover, the evolution of their physical properties with 
increasing cluster size provides new perspectives for un- 
derstanding the microscopic origin of condensed-matter 
properties [j]]. In this context, the study of magnetism 
in metals posses a particularly interesting and challeng- 
ing problem since the properties of magnetic metals of- 
ten change dramatically as the electrons of an isolated 
atom become part of a cluster of several atoms and de- 
localize. It is therefore of fundamental importance to 
understand how itinerant magnetism, as found for exam- 
ple in 3d transition-metal (TM) solids, develops starting 
from localized atomic magnetic moments. The strong de- 
pendence of the magnetic behavior as a function of sys- 
tem size, structure and composition opens in addition the 
possibility of using clusters to tailor new magnetic ma- 
terials for specific technological purposes. Consequently, 
the relation between local atomic environment, cluster 
structure and magnetism has also implications of practi- 
cal interest. 

Most theoretical studies of itinerant cluster magnetism 
have been performed using Hartree-Fock (HF) and local 
spin density (LSD) methods g. Exact many-body cal- 
culations are presently limited to simple models, such as 
the Hubbard model, and to systems containing a small 
number of sites |}|-|| . Within mean-field approximations 
(HF or LSD) the self-consistent spin polarizations are 



usually restricted to be collinear, i.e., the direction of all 
local magnetic moments is assumed to be the same. How- 
ever, the well known sensitivity of itinerant magnetism 
to the local environment of the atoms Q suggests that 
other instabilities towards spiral-like spin-density-waves 
(SDW's) (|] or even more complex magnetic structures 
should also be possible in general. For example, the ten- 
dency to antiferromagnetic (AF) order close to half-band 
filling and the presence of triangular loops in compact 
structures generates magnetic frustrations that may eas- 
ily yield complex arrangements of the local magnetic mo- 
ments. Moreover, the reduction of local coordination 
numbers at the surface of clusters removes constraints 
between the local moments and could favor the devel- 
opment of low-symmetry spin polarizations. In solids, 
noncollinear magnetic structures have been identified ex- 
perimentally and theoretically already for a long time 
In contrast, very little is known at present in the 
case of finite systems pO|-p^| . Fully unrestricted ab initio 
calculations of noncollinear spin arrangements have been 
performed only recently for very small Fe^ clusters [ p"2| . 
The investigation of magnetic phenomena of this kind 
requires a symmetry unrestricted approach in which no 
a priori assumptions are made concerning the relative 
orientation of the local magnetic moments, thereby en- 
larging the number of degrees of freedom of the problem. 

The main purpose of this paper is to investigate the 
characteristics of noncollinear magnetic states in finite 
clusters. We consider the single-band Hubbard model 
and determine the ground-state magnetic properties in 
a fully unrestricted Hartree-Fock approximation. The 
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theoretical approach, outlined in Sec. [n[ is applied to 
clusters having N < 43 atoms. The results presented 
in Sec. Ill analyze noncollinear magnetic behaviors for 
a few representative compact structures. Several ex- 
amples are given that illustrate the large variety of 3- 
dimensional magnetic arrangements obtained in the self- 
consistent calculations as a function of the Coulomb re- 
pulsion strength U/t, band filling v/N and cluster struc- 
ture. In Sec. |fv| we compare the UHF results with exact 
diagonalization calculations for N < 13 atoms (Lanczos 
method) j5j. The role of quantum fluctuations beyond 
mean field and the consequences of the often artificial 
breaking of spin-symmetry implied by the formation of 
the noncollinear local moments are discussed. Taking 
into account that approximations such as UHF are un- 
avoidable for larger clusters and for more realistic model 
Hamiltonians, it is of considerable interest to test the 
validity of these methods in order to improve the inter- 
pretation of the approximate results and to obtain a more 
accurate description of magnetic phenomena in clusters. 



II. THEORETICAL METHOD 



The single-band Hubbard Hamiltonian |l3| is given by 



H 



(l,m) ,er 



U J2 



(2.1) 



The single-particle states 



(2.3) 



are linear combinations of the atomic-like orbitals asso- 
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Notice that in Eq. fl2.3|) we allow for the 



most general superposition of single-electron states since 
al may involve a mixture of both | and J, spin com- 
ponents. The eventually complex coefficients A\ a are 
determined by minimizing the energy expectation value 
£"uhf = (UHF|_ff |UHF). In terms of the density matrix 



PUTtma , = (UHF|4 w |UHF) = AtA k ma , , (2.4) 



k=l 



this is given by 



-E"uhf 



(2.5) 



The energy minimization and the normalization con- 
straints on the wave function lead to the usual self- 
consistent equations 



- t 



; A k ia + U {A k laPlsra - A k sP i sM ) = e k A k a . (2.6) 



where c[ a (q<j) is the creation (annihilation) operator of 
an electron at site I with spin a, and hia = c\ a ci a is 
the corresponding number operator [fii — X^o-'W)- The 
first sum runs over all pairs of nearest neighbors (NN) 
and the second over all sites. The model is character- 
ized by the dimensionlcss parameter U/t, that measures 
the relative importance between kinetic and Coulomb en- 
ergies, by the cluster structure, that defines the kinetic 
energy operator, and by the number of valence electrons 
v. The variations of U/t can be associated to a uniform 
relaxation of the interatomic distances (e.g., t oc R^ 5 
for TM's) or to changes in the spatial extension of the 
atomic-like wave function, as in different elements within 
the same group. Different z/s correspond to different 
band-fillings v/N that may be associated qualitatively 
to the variations of v/N across a TM d series. In spite of 
its simplicity, this Hamiltonian has played, together with 
related models, a major role in guiding our understand- 
ing of the many-body properties of metals and of low- 
dimensional magnetism. It is the purpose of this work to 
use it to investigate the properties of noncollinear itiner- 
ant magnetism in small compact clusters. 

In the unrestricted Hartree-Fock (UHF) approxima- 
tion the ground state for v electrons is a single Slater 
determinant that can be written as 



|UHF) 



n« 



L k 



(2.2) 



For a given solution, the average local electronic density 
(ni) is given by 



(«;} = Pi],i] + Pn,n > 

and the spin polarizations vectors (Si) 
({Sf),(Sf),{Sf)) by 

(S?) = ( Pim + Pll<n )/2, 
(Sj) = -i (pi Ul -pi Ur )/2, 
(S?)= (p nA -p im )/2. 



(2.7) 



(2.8) 



The usual collinear UHF approach is recovered when 
Pio,i5 = V I, i.e., when all magnetic moments (Si) are 
parallel to z. In practice, several random spin arrange- 
ments are used as starting points of the selfconsistent 
procedure in order to ensure that the final result corre- 
sponds to the true UHF ground state. In case of multiple 
selfconsistent solutions (nonequivalent by rotations) the 
UHF energies are compared. -Euhf can be rewritten as 

£ UHF = -t Y, Pl*,m* + ^Y,(ni) 2 -Uj2\(Si)\ 2 . 



2_ Pl<r,ma + J 2^ n '/ 
(l,m),cr 



(2.9) 



One observes that the Hartree-Fock Coulomb en ergy 
the sum of the 2nd and 3rd terms in Eq. (2.9) 
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— favors a uniform density distribution and the forma- 
tion of local moments (Si). Due to the local character of 
Hubbard's Coulomb interaction, the relative orientation 
of the different (Si) does not affect E^ F . It is therefore 
the optimization of the kinetic energy that eventually 
leads to the formation of complex magnetic structures 
with | (Si) ■ (S m ) | ^ 1 or to non-uniform density distribu- 
tions (n{). 

As a result of the tendency to avoid double orbital oc- 
cupancies, the UHF solutions often correspond to states 
of broken symmetry: spin-density waves (SDW's), charge 
density waves (C DW ' s) o r both. The spin-rotational in- 



variancc of Eqs. (2.2-2.6) implies that the energy is un- 
changed after a rotation of the whole spin arrangement 
{(Si), I = 1,...,N}. Therefore, if (Si) + one has a 
set of linearly independent congruent solutions UHFfc) 
(k > 2) which have the same average energy -Euhf and 
which differ from each other only by the orientation of 
the spin polarizations relative to the cluster structure. 
The illustrations of spin arrangements shown in Sec. Ill 
correspond to one of these SDW's, which is chosen only 
for the sake of clarity. The set of UHF solutions may 
be used to restore the symmetry appropriate to the ex- 
act ground state \^o) thereby improving the approxi- 
mate wave function. For instance, for some SDW's hav- 
ing Y)r(Si) — one may consider the spin-symmetrized 
Hartree-Fock (SSHF) wave function |iHJT~ 



|SSHF) = 



|UHF1) + |UHF2) 
v/2(l + (UHF1|UHF2)) 



(2.10) 



where |UHF2) is obtained form |UHF1) by interchang- 
ing up and down spins. In spite of its simplicity, |SSHF) 
goes beyond the UHF approximation and corresponds to 
a correlated state satisfying 



|(SSHF|* ) 



= ^2/(1 + (UHF1|UHF2)) |(UHF1|# )| 
> |(UHF1|* )| - |(UHF2|tf >| • (2-11) 



If quantum fluctuations between the two SDW's |UHF1) 
and |UHF2) are non-negligible ((UHFl|i? |UHF2) ^ 
(UHF1|UHF2)_Euhf) an energy reduction relative to 
i?UHF is obtained since 



E' 



SSHF 



£uHF + (UHFl|iI|UHF2) 
1 + (UHF1|UHF2) 



(2.12) 



AE = Euhf — -Esshf measures the importance of these 
quantum spin fluctuations. More complex symmetrized 
states involving a linear combination of several degener- 
ate UHF states may be constructed analogously. 

In order to quantify the role of electron correlations 
and to assess goal and limitations of the UHF approx- 
imation in applications to finite clusters, we shall also 
compare some of our results with those obtained by ap- 
plying exact diagonalization methods (^-|^| . In this case 
the Hubbard model is solved numerically by expanding 
its ground-state |*o) m a complete set of basis states 



|$ m ) which have definite occupation numbers n^J: at all 
orbitals la (ri iCT |$ m ) = <|$ m ), with nf a = 0, 1). |tf„) 
is written as 



where 



l$m) 



n 

la 



(4)"" 



\vac) 



(2.13) 



(2.14) 



The values of satisfy the usual conservation of the 
number of electrons v = + v j and of the z component 
of the total spin S z — (i/j — v{)/2, where v a — J^i n uy 
Taking into account all possible electronic configurations 
may imply a considerable numerical effort which depends 
on the number of atoms and on band filling. For not 
too large clusters, the expansion coefficients ai m can be 
determined by sparse-matrix diagonalization procedures. 
The results presented in this work were obtained using a 
Lanczos iterative method j|,|]jl(|. In order to calculate 
|$o) one usually works in the subspace of minimal S z 
since this ensures that there are no a priori restrictions 
on the total spin S. The ground-state spin S is then 
obtained by applying to |\&o) t ne total spin operator 



s 2 — ^ §i ■ s m — 

Irn lm L 



(2.15) 



where S^~ 



c] t qj., S t = c^Q-f and Sf = (njf - "u)/2. 

From the expectation values of Si ■ S m one also obtains 
the local magnetic moments (I — m) and the intersite 
spin correlation functions (I ^ m). 



III. NONCOLLINEAR MAGNETIC ORDER IN 
COMPACT CLUSTERS 

In Fig. | UHF results are given for the local magnetic 
moments /if = \(Si)\ and the total magnetic moment 

Mt = I Y^i(@l)\ °f fcc-like clusters having N — 13, 19 and 
43 atoms at half band filling (y — N) . These clusters are 
formed by adding to a central atom (/ = 1) the successive 
shells of its first NN's (N = 13), second NN's (N = 19) 
and third NN's (N = 43). Several common properties 
are observed as a function of U jt. Starting from the un- 
corrected limit (U = 0), the total moment /it remains 
approximately constant for U/t < 3-4. In the weakly 
interacting regime, the local moments \xi do not depend 
strongly on U/t and are generally small. Notice that for 
N = 43, fir is not minimal at small U (/Xt = 3/2) due to 
degeneracies in the single-particle spectrum. For larger 
U/t, \it decreases rapidly eventually with some discon- 
tinuities and a few oscillations reaching values close to 
Ht = for U/t ~ 5-6 (U/t ~ 9 for N — 43). At this 



3 



intermediate range of U/t, the \ii increase more rapidly 
reaching values not far from saturation at the U/t for 
which fix is minimum (//; ~ 0.40-0.45). The opposite 
trends shown by fii and \xt are a clear indication of the 
expected onset of strong AF-like order at half-band fill- 
ing. We shall see in the following that this corresponds 
in fact to noncollincar spin arrangements. If U/t is in- 
creased beyond U/t = 5-6 (beyond U/t — 9 for N = 43) 
the local moments do not vary significantly, and [It ei- 
ther does not change very much (N = 19) or increases 
monotonically (N = 13 and 43) remaining always smaller 
than 1/2. 

The changes in fii and fix are the result of qualitative 
changes in the magnetic order. As an example we show 
in Fig. ^| the selfconsistent spin arrangements obtained in 
fee clusters with N = 13 atoms for representative values 
of U/t (v = N) @- For small U/t (U/t < 3.7) one finds 
a collinear AF order with small \ii [Fig. ||(a)]. Here we 
observe a charge and spin density-wave at the cluster sur- 
face that is related to degeneracies in the single-particle 
spectrum at U = 0. The atoms belonging to the central 
(001) plane [shown in grey in Fig. ||(a)] have much larger 
moments than the atoms at the upper and lower (001) 
planes. For example for U/t = 0.5, /i; = 0.26/xb at the 
the central plane, while the other surface moments are 
/if = 0.0046/xb (see Fig. |l|). The central atom has a small 
spin polarization (hi — 0.002/^b for U/t — 0.5). No- 
tice that the magnetic order within the upper and lower 
(001) planes is ferromagnetic- like (with small fii) and 
that the surface moments belonging to successive (001) 
planes couple antiferromagnetically. Thus, the magnetic 
moments at the surface of the central plane are not frus- 
trated since all its NN's are antiparallel to them. This 
explains qualitatively the larger /i; found at this plane. 
In contrast, unavoidable frustrations are found for the 
smallest magnetic moments at the central site and at the 
atoms of the upper and lower (001) planes [see Fig. ||(a)] 

The crossover from the smal\-U/t to the large-?7/£ 
regime takes place as a succession of noncollinear spin 
arrangements which attempt to minimize the magnetic 
frustrations among the increasing local moments. A rep- 
resentative example is shown in Fig. |^(b). While in this 
case the spin arrangement is noncollinear, all the spin 
moments still lie in the same plane. The very small val- 
ues of [It shown in Fig. [I] for intermediate U ft indicate 
that there is an almost complete cancellation among the 
(Si). However, this type of spin arrangements are quite 
unstable if the strength of Coulomb interactions is fur- 
ther increased. At U/t ~ 5.1 the cluster adopts a fully 3- 
dimensional spin structure which remains essentially un- 
changed even in the strongly correlated limit [Fig. @(c)]. 
The spin arrangement can be viewed as a slight distor- 
tion of the spin ordering that minimizes the energy of a 
classical AF Heisenberg model on the surface shell, i.e., 
ignoring the interactions with the central site. In fact, 
if the central atom were removed or if it carried no lo- 
cal moment, as it is the case for v — 12, the surface 



moments (Si) would point along the medians of one of 
the triangles at the surface and would lie all within a 
plane. The magnetic interactions with the central spin 
(S\) in the 13-atom cluster induce a small tilt (Sf) of the 
surface spin polarizations (Si) which is opposite to (Si) 
(0l) || z for v = N = 13). The (Sf) component is the 
same for all surface sites and depends moderately on U/t 
(\(Sf)\ = 0.023-0.056 for U/t > 5.1). 

Similar magnetic structures are also found in larger 
symmetric fee clusters. In Fig. || the self-consistent spin 
configuration for N — 43 and U/t — 10 is illustrated 
(v = N). The moment (Si) of the central site points 
along the (111) direction (± to the plane of Fig. |) [|7). 
The other (S m ) lie almost in the plane of the figure 
with only small components (S^) along the (111) direc- 
tion. As for N = 13, (S^j) is induced by the interac- 
tions with (Si). In fact, if (Si) vanished, all the spins 
would be in the plane of the figure. (S^) is the same 
for all atoms in a given NN shell, and changes sign as 
we move from the center to the surface of the cluster. 
For example for U/t = 10, cos6q TO = —0.22 for m be- 
longing to the first shell, cos#i m = 0.25 for m in the 
second shell and cos#i m = 0.01 for m in the third shell 
(cos6» im = (Si) ■ (S m )/\(Si)\\(S m )\). Notice that the spin 
ordering of the innermost 13 atoms is similar to the one 
found in fee clusters with N = 13 [Fig. §(c)]. While 
these trends hold for symmetric fee clusters, it is worth 
to remark that strong modifications of the magnetic order 
generally occur if the symmetry of the cluster is lowered. 
For instance, for N = 14 — an atom added to the closed- 
shell 13-atom cluster — one obtains a ferromagnetic-like 
coupling between the central atom and some of the first 
neighbors (v = N and U/t = 10). 

A more detailed account of the spin correlations in the 
UHF ground state of fee clusters is provided by the ex- 
pectation values (Si-S m ). In Fig. ||the average of (Si-S m ) 
between NN atoms at different shells is shown. In most 
cases we observe that (Si-S m ) < (AF correlations) and 
that |(5z-5 m )| increases with increasing U/t. A particu- 
larly important increase of AF correlations is observed for 
intermediate U/t when the local moments fii are formed 
(U/t ~ 4.5-5.5 for N = 13 and 19, see also Fig. Q). This 
is consistent with the already discussed decrease of fix 
and the onset of AF order. There are, however, two ex- 
ceptions to this trend. The correlations 701 between the 
central site and the surface shell for N — 13 increase 
somewhat for U/t ~ 5 (701 < 0) implying that these AF 
correlations first tend to be less important when the /i/ 
increase. It seems that, as a result of frustrations, the 
more important increase of AF correlations 711 among 
the surface NN's is done at the expense of the correla- 
tions between the central atom and the surface. A similar 
behavior is found in exact calculations, as it will be shown 
in Sec. IV. Another interesting case concerns the corre- 
lations 712 between the first and second atomic shells for 
N = 19 and 43. Here we observe that 712 changes sign 
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upon going from N = 19 to N — 43 (U/t > 5). More- 
over, once the local moments are formed, 701 decreases 
for N = 19 [701(19) < 0] and increases for N — 43 as 
U/t increases [701 (43) > 0]. This behavior can be qual- 
itatively understood by comparing the surfaces of these 
clusters. For N = 19 the outermost spins (shell 2) are 
free to couple antiferromagnetically with their NN's (shell 
1). However, the presence of an additional atomic shell 
for N = 43, which has NN's belonging to both the first 
and second shells, forces a parallel alignment of the spins 
in shells 1 and 2 in order to allow AF coupling with the 
third shell. In fact, as shown in Fig. ||(c), 7i3 an d 723 
are the strongest AF correlations in the 43-atom cluster. 
The same conclusions are drawn by comparing the rel- 
ative orientations of the spin polarizations (Fig. |^). In 
particular we observe that the angles 823 between NN 
spin polarizations at shells 2 and 3 are the largest in av- 
erage (623 — 7r or close to it). 

In addition to the AF half-filled case it is also interest- 
ing to discuss other band fillings, for example v = N + 1 
that is known to develop a FM ground state in the limit 
of large U/t @. In Fig. | results are given for \xi and 
fj,T in an fee 13-atom cluster with v = 14 electrons. In 
this case we observe an essentially monotonic increase of 
Ht that reflects the progressive development of a fully 
polarized ferromagnetic state. Close to the threshold 
for the onset of ferromagnetism (4.5 < U/t < 6.5) the 
changes of the AF-like spin arrangement produce small 
oscillations of fix (see Fig ||). The approximately step- 
like behavior resembles at first sight the results obtained 
in collinear mean-field calculations p9[ . However, in the 
present case the physical picture behind the formation of 
a FM state is quite different. The increase of [It close 
to the steps involves a succession of noncollinear spin 
arrangements with increasing degree of parallel moment 
alignment. Moreover, notice that the local moments in- 
crease much more rapidly than \ix approaching satura- 
tion (hi ~ 0.40-0.45) for values of U/t at which \ix is still 
small. Thus, the increase of fix with U/t is the result of 
the parallel alignment of already existing local moments 
fil- In contrast, in collinear Hartree-Fock calculations 
the increase of /it is associated to the formation of the 
local moments themselves since fix approximately pro- 
portional to [ii p9[ | . Comparison with exact calculations 
shows that the present noncollinear picture is qualita- 
tively closer to the actual ground-state magnetic prop- 
erties than the one derived in collinear calculations, at 
least for the single-band Hubbard model. Indeed, we 
have computed the local moments nj — (Sf) and the 
total spin S of the same fee 13-atom cluster (v = 14) 
as a function of U/t by using a Lanczos exact diago- 
nalization method. Starting from the uncorrelated limit 
[f4(U = 0) = 0.28 and fif(U = 0) = 0.36 for I = 2 
13)] fxf increases rapidly with U/t reaching values close 
to saturation already for U/t ~ 10 \pl(U = 10) = 0.56 
and fj,f(U = 10) = 0.64 while n\{U = 00) = 0.56 and 
fjtf(U = 00) = 0.70 (I = 2-13)]. In contrast the ground- 
state spin S remains equal to zero up to U/t ~ 40. This 



implies that for v — N + 1 = 14 the local moments are 
formed well before FM order sets in, as observed in the 
noncollinear UHF calculations (Fig. ||). Notice, however, 
that the situation could be different for other band fill- 
ings where ground-state ferromagnetism is found at much 
smaller U/t or where S is a non-monotonous function of 
U/t (e.g., for v = 15-18 f§). UHF yields larger local 
magnetic moments at the cluster surface in agreement 
with the exact results but it underestimates severely both 
the Coulomb repulsion U/t above which /xt > 1/2 and 
the U/t for reaching saturation. 

Calculations have been also performed for fcc-like 13- 
atom clusters having ^ = ^ — 1 = 12 electrons. For small 
U (U/t < 4.9) the calculated magnetic order is collinear 
((Si) || z). The UHF ground state is a broken symmetry 
state that shows a charge- and spin-density wave along 
the (001) direction. For U/t < 3.5 the atoms I in the 
central (001) plane present an enhanced electron density 
(ni) ~ 1.12 and local magnetic moments /i; = 0. The 
atoms of the upper and lower (001) planes show AF or- 
der within each plane and \ii = 0.17. The vanishing fii 
at the central (001) plane can be interpreted as a con- 
sequence of magnetic frustrations since the sum of their 
NN spins is zero. For U/t > 3.5, the surface atoms at 
the central (001) plane develop local moments \ii 7^ 0. 
The spin arrangement remains collinear with important 
frustrations at some triangular faces. For U/t > 4.9 the 
spin-density distribution changes to noncollinear with all 
surface spin having the same modulus. However, the cen- 
tral site remains unpolarized. In fact the magnetic order 
is the same as the one obtained if the central site is re- 
moved (N — 12). The moments at the triangular faces 
point along the medians just as in the classical antifer- 
romagnetic ground state of an isolated triangle (TTJ . In 
spite of the three-dimensional geometry of the cluster the 
spin structure is two-dimensional (all magnetic moments 
lie on the same plane). 

The topology and symmetry of the cluster structure 
plays a major role in determining the magnetic order and 
magnetic correlations within the cluster. In this context, 
it is interesting to consider different geometries and to 
compare their magnetic behavior. In Fig. |] the UHF 
magnetic order in a 13-atom icosahedral cluster having 
v = 12 electrons is illustrated for representative values 
of U/t. As a result of degeneracies in the single-particle 
spectrum (U — 0) a noncollinear, coplanar arrangement 
of the local magnetic moments is found at the surface 
already for very small U/t (U/t < 5, see Fig. ^). At 
the central atom (not shown in the figure) the magnetic 
moment fi% = 0. A similar behavior is observed in fee 
clusters (N = 13, v — 12). The tendency to avoid frus- 
trations within NN triangular rings conditions the ar- 
rangement of the local moments. In some triangles the 
spin polarizations point along the medians just as in the 
classical Heisenberg model, while in others two spins are 
parallel (i.e., fully frustrated). For larger values of U/t 
(U/t > 5) a less frustrated solution is favored by the an- 
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tiferromagnetic correlations. This corresponds to a truly 
three-dimensional arrangement of the surface moments 
with hi — 0. The spin structure is such that if the sur- 
face (Si) are brought to a common origin they form an 
icosahedron. The magnetic moments (Si) on pentagonal 
rings present a small component that is perpendicular to 
the plane containing the atoms and that is antiparallel 
to the magnetic moment of the atom capping the ring. 
The projections of (Si) on to the plane of the ring are 
ordered in the same way as in an isolated pentagon. 



IV. COMPARISON BETWEEN UHF AND 
EXACT RESULTS 

In Fig. [?] UHF and exact results are given for the lo- 
cal moments /if = (Sf) and spin correlations (Si ■ S m ) 
in an fcc-like 13-atom cluster. The UHF results for /sf 
are quantitatively not far from the exact results. Not 
only the uncorrelated limit (U = 0) is reproduced, but 
good agreement is also obtained in the large U / t regime. 
Main trends such as the larger jif at the cluster surface 
for small U/t and the reduction of the difference between 
surface and inner moments for large U/t are correctly 
given. However, UHF underestimates the increase of [ij 
for U/t < 3.7 and anticipates the tendency to localization 
with increasing U/t. This results in a much more rapid 
crossover from weak to strong interacting regimes than 
in the exact solution. The quantitative differences are 
more important in the case of spin correlation functions 
(Si-Sm) particularly for large U/t [Fig. 0(b)]. Here we 
find that UHF underestimates the strength of AF spin 
correlations 711 at the surface, since the formation of 
permanent local moments blocks quantum spin fluctua- 
tions along the transversal directions. Still, in both UHF 
and exact calculations, the increase of AF correlations at 
the surface (increase of |7n|) is done at the expense of 
a decrease of the spin correlations with the central atom 
(decrease of |7oi|)- The discrepancies between UHF and 
exact results for N = 13 show the limits of mean-field 
and give us an approximate idea of the corrections to be 
expected in correlated calculations on larger clusters. As 
expected, UHF yields better results for properties like 
the local moments, that are related to the density distri- 
bution, than for the correlation functions. Nevertheless, 
since the trends given by UHF are qualitatively correct, 
it is reasonable to expect that the conclusions derived for 
larger clusters are also valid (Figs, [l] and ^). 

The determination of the structure of magnetic clusters 
is a problem of considerable importance since structure 
and magnetic behavior are interrelated Moreover, 
since most calculations of cluster structures are based 
upon mean field approximations, it would be very inter- 
esting to evaluate the role of electron correlations. We 
have therefore determined the relative stability of a few 
representative cluster structures as a function of band 
filling v/N and Coulomb repulsion strength U/t in the 



framework of the UHF approximation to the Hubbard 
model and we have compared the results with available 
exact calculations || . Four different symmetries are con- 
sidered: icosahcdral clusters, which maximize the aver- 
age coordination number, face centered cubic (fee) and 
hexagonal close packed (hep) clusters, as examples of 
compact structures which are found in the solid state, 
and body centered cubic (bec) clusters, as an example of 
a rather open bipartite structure. These cluster geome- 
tries are representative of the various types of structures 
found to be the most stable in rigorous geometry opti- 
mizations for A < 8 The results for N = 13 are 
summarized in the form of a magnetic and structural 
diagram shown in Fig. ||. A qualitative description of 
the type of magnetic order obtained in the self-consistent 
calculations is indicated by the different shadings. One 
may distinguish three different collinear spin arrange- 
ments: non-magnetic solutions (NM), non-saturated or 
weak ferromagnetic solutions (WFM) and saturated fer- 
romagnetic solutions (SFM). The NM case includes para- 
magnetic states in which the total moment \iT is minimal 
(Ht = or 1/2). Concerning the noncollinear spin ar- 
rangements we distinguish two cases: noncollinear non- 
magnetic states (NC) in which non-vanishing (eventu- 
ally large) local moments fii sum up to an approximately 
minimal total moment fiT < 1, and noncollinear ferro- 
magnetic states (NCFM) that show a net magnetization 
Mr > 1- The NC states include all sort of frustrated 
antiferromagnetic-like spin structures, for example, those 
illustrated in Figs. || and ^. In order to quantify the effect 
of electron correlations we also show in Fig. ^| the corre- 
sponding magnetic and structural diagram as recently 
obtained by using exact diagonalization methods || . 

For small U/t (U/t < 10) the UHF results for the most 
stable of the considered structures are in very good agree- 
ment with the exact calculations [compare Figs. || and |). 
The icosahedral cluster yields the lowest energy in the low 
carrier-concentration regime (v < 6), as could have been 
expected taking into account that the largest coordina- 
tion number is favored for small v as in smaller clusters 
Q. In this case, the kinetic energy dominates (uncorre- 
lated limit) and therefore the structure with the largest 
bandwidth is stabilized. As v is increased several struc- 
tural transitions occur. For small U/t, both exact and 
UHF calculations present the same structural changes: 
from icosahedral to fee structure at v = 11, from fee to 
hep at v = 17, and from hep to bec at v = 20. At 
larger U/t (U/t > 12) the interplay between the kinetic 
and Coulomb energies introduces important correlations 
that cannot be accounted for within UHF and that play 
a central role in the determination of the magnetic and 
structural properties. Thus, UHF often fails to yield the 
lowest-energy structure in the limit of large U/t. A main 
source of discrepancy is the too strong tendency of UHF 
to yield SFM ground states, particularly above half-band 
filling, which often disagrees with the exact magnetic be- 
havior. Consequently the optimal structure is missed 
rather frequently. For example, for v — 19 and large 
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U/t, UHF predicts the fee structure with a SFM ground 
state, while in the exact calculation the icosahedral clus- 
ter with S = 1/2 is the optimum. Similar drawbacks are 
seen for other band fillings such as v = 7, 10, 21, and 22 
(large U/t). Still, in the event that the true ground state 
does show strong ferromagnetism, UHF succeeds since 
the ground state is the superposition of single-particle 
states. Examples of this kind are v = 12 and v — 20 for 
large U/t. These are rather exceptions, however, since 
the presence of a FM ground state in the exact solution 
is in general much less frequent than predicted by UHF. 
Around half-band filling UHF reproduces qualitatively 
well the transition from fee or hep to icosahedral struc- 
ture with increasing U/t as well as nontrivial re-entrant 
effects (y = 13-18). However, UHF underestimates the 
value of U/t at which the structural changes occur (see 
Figs. H and ||) . Part of the quantitative differences could 
be remedied by using in the UHF calculations a reduced 
or renormalized U/t which simulates some of the effects 
of correlations. 
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Summarizing, UHF fails to reproduce the exact phase 
diagram in detail, particularly in some of the most inter- 
esting antiferromagnetic or weak ferromagnetic regimes. 
Structural transition are sometimes missing (e.g., for 
v = 10) and in other cases changes of structures appear 
artificially (e.g., bec to hep for v — 21 and 22). Nev- 
ertheless, it is also fair to say that UHF yields a good 
account of the relative stability between the considered 
structures well beyond the weakly interacting limit (up 
to U/t ~ 16) and that it also explains the larger stability 
of ferromagnetism above half-band filling. In the limit of 
strong interactions the validity of UHF breaks down and 
correlation effects beyond the single-determinant wave- 
function are indispensable in order to obtain the ground- 
state structure and magnetic behavior reliably. Improve- 
ments on the UHF wave function could be introduced by 
restoring the broken spin symmetry as proposed in Refs. 
QUI [Eqs. ( |2.10| - |2.12D ]. However, the success of such an 



approach is likely to depend on the geometry of the clus- 
ter. For example, in low symmetry structures UHF tends 
to exaggerate the formation of spin and charge density 
redistributions which may be far from the actual exact 
solution and which can be very difficult to restore a pos- 
teriori. In such cases a Jastrow-like variational ansatz 
on restricted Hartree-Fock states could be more appro- 
priate. Finally, it should be recalled that the Hubbard 
model for small clusters with one orbital per site is proba- 
bly the most extremely low-dimensional system one may 
consider. Therefore, fluctuations and correlations effects 
are expected to be here more drastic than in larger clus- 
ters or in realistic multiband models more appropriate 
for the description of TM's. 
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FIG. 1. Magnetic moments as a function of the Coulomb 
repulsion strength U/t for fee-like clusters having (a) N = 13, 
(b) N = 19 and (c) N = 43 atoms at half-band filling (UHF 
approximation). Local moments = |(Si)| are shown for the 
central atom I = 1 (dashed), its first NN's (I = 2-13, dot- 
ted), second NN's (/ = 14-19, dashed-dotted) and third NN's 
(I = 20-43, short dashed). In (b) and (c) refers to the shell 
average. The total moment \xt = 1X^(^)1 ^ s gi ven by the 
solid curves. 
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FIG. 3. Illustration of the local magnetic moments (Si) 
in a 43-atom fee- like cluster at half-band filling and U/t = 10 
As in Fig. |(c), the plane of the fi gure is perpendicular 
to the (111) direction. The atoms are represented by circles 
using different sizes and grey tones for different (111) lay- 
ers. The later are denoted by A, B, C and A' starting from 
above. The arrows show the projection of (Si) on to the (111) 
plane. At the central atom {Si) || (111). The spin structure 
is 3-dimensional with off-plane components of (Si) that alter- 
nate sign on different shells around the central atom. For the 
13 innermost atoms the magnetic order is similar to the one 
shown in Fig. |(c). 



FIG. 2. Illustration of the UHF magnetic order in fee-like 
clusters having N = 13 atoms, v — N — 13 electrons and rep- 
resentative values of U/t (Hubbard model) @. (n) U/t < 3.7: 
The surface atoms in grey have larger local magnetic moments 
(5;) | = 0.16-0.18 and less electron density {m) ~ 0.88-0.90 
than the other surface atoms for which \{Si)\ — 0.02-0.03 
and (ni) — 1.05-1.12. The spin arrangement is collinear. (b) 
4.8 < U/t < 5.1: The plane of the figure is perpendicular to 
the (111) direction. The arrows show the local spin polariza- 
tions {Si) which are all within the (111) plane, (c) U/t > 5.1: 
The plane of the figure is perpendicular to the (111) direction. 
The arrows show the projection of (Si) on to the (111) plane. 
At the central atom (Si) || (111). The off-plane components 
of (Si) at the surface atoms are antiparallel to the central 
spin and have all the same value (|(Sf}| ~ 0.023-0.056 for 
I = 2-13). 
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FIG. 4. Spin correlations {Si ■ S m ) in fcc-like clusters hav- 
ing N = 13-43 atoms as a function of U jt. Results are given 
for the average spin correlations 701 between the central site 
and its first NN's (solid), 711 between NN sites on the first 
shell (dashed), 712 between NN sites on the first and second 
shell (dashed-dotted) , 713 between NN sites on the first and 
third shell (dotted), and 723 between NN sites on the second 
and third shell [lower thin curve in (c)]. 



FIG. 5. Magnetic moments in an fcc-like 13-atom cluster 
with v — 14 electrons as a function of U/t (UHF approxima- 
tion): (a) total moment /it = | anc3 - 3 ) l° ca l moments 
fJ-i = I (Si) I at the central site (dashed line) and at the surface 
(dotted line). 
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FIG. 7. Comparison between UHF and exact results for 
an fcc-like 13-atom cluster at half-band filling, (a) Local mag- 
netic moments fif — \(SrSi)\ at the central site (1 = 1, dashed 
curves) and at the cluster surface (/ = 2-13, solid curves), 
(b) Average 7^ of the spin correlation functions {SrS m ) be- 
tween the central site and its first NN's at the surface (701, 
dashed curves) and between NN's at the surface shell (711, 
solid curves). The curves with (without) crosses correspond 
to exact (UHF) results. 



FIG. 6. Illustration of the magnetic order at the surface 
of an icosahedral 13-atom cluster with v = N — 1 = 12 elec- 
trons (UHF approximation) @. (a) 0.8 < U/t < 5 : All 
local magnetic moments (Si), represented by the arrows, are 
parallel to the plane of the figure, (b) U/t > 5 : The arrows 
show the projections of {Si} on the plane of the figure. At site 
1 — 2, (Si) is perpendicular to the plane of the figure (along 
the z axis). The (Si) at sites I = 3-7 have antiparallel projec- 
tion to (S 2 ) [cose 2l = ((St) ■ (&»/(|(#)|K&)| ~ -CU5)]. If 
brought to a common origin, the spin polarizations (Si) form 
a perfect icosahedron. Notice that the positions of the atoms 
in the lower pentagonal rings have been expanded along the 
polar radius in order to ease the visualization. 
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FIG. 8. Magnetic and structural diagram of clusters 
having N — 13 atoms as obtained by using the Hubbard 
model in the unrestricted Hartree-Fock (UHF) approxima- 
tion. Four different types of structures are considered: icosa- 
hedral (ICOS), face centered cubic (FCC), hexagonal close 
packed (HCP), and body centered cubic (BCC). The low- 
est-energy structure is given as a function of the Coulomb 
repulsion strength U, hopping integral t and number of elec- 
trons v. The bold lines indicate the structural transitions. For 
v < 6 the lowest-energy structure is icosahedral and for v > 23 
it is bcc. The corresponding UHF magnetic orders are indi- 
cated by different shading patterns: saturated ferromagnetic 
(crossed lines, SFM), collinear weak ferromagnetic (diagonal 
lines, WFM), collinear nonmagnetic (no shading, NM), non- 
collinear with total moment fir < 1 (horizontal lines, NCT), 
and noncollinear with total moment \it > 1 (vertical lines, 
NCFM). 
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FIG. 9. Magnetic and structural diagram of Hubbard clus- 
ters having TV = 13 atoms as obtained by using exact diago- 
nalization methods (after Ref. ||). As in Fig. g, icosahedral 
(ICOS), face centered cubic (FCC), hexagonal close packed 
(HCP) and body centered cubic (BCC) structures are con- 
sidered. The most stable structure is given as a function of 
Coulomb repulsion U, hopping integral t and number of elec- 
trons v. The corresponding ground-state spin S is minimal 
(S = or 1/2) unless explicitly indicated. Broken lines sep- 
arate regions having the same structure but different S. For 
v < 9 the icosahedron yields the lowest energy for all U/t, 
while for v > 21 the bcc structure is the most stable. 



13 



